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Using femtosecond pump-probe spectroscopy, we identify excitation induced dephasing as a major 
mechanism responsible for the breaking of the strong-coupling between excitons and photons in a 
semiconductor microcavity. The effects of dephasing are observed on the transmitted probe pulse 
spectrum as a density dependent broadening of the exciton-polariton resonances and the emergence 
of a third resonance at high excitation density. A striking asymmetry in the energy shift between 
the upper and the lower polaritons is also evidenced. Using the excitonic Bloch equations, we quan¬ 
tify the respective contributions to the energy shift of many-body effects associated with Fermion 
exchange and photon assisted exchange processes and the contribution to collisional broadening. 
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A semiconductor microcavity is a system that confines 
photons and allows them to strongly interact with quan¬ 
tum well excitons [ 1 ]. Polaritons are composite parti¬ 
cles arising from the coherent superposition of a pho¬ 
ton and exciton. The interactions mediated by their 
excitonic part make semiconductor microcavities a suit¬ 
able platform for realizing nonlinear optical devices such 
as: bistablity memory [2] and polariton switching [3]. 
Moreover, applications of semiconductor microcavities to 
quantum information have been also proposed on the ba¬ 
sis of their coherent behaviour H E]. When modeling 
these devices, they are usually described with a nonlin¬ 
ear Schrodinger equation, which is formally the same as 
the Gross-Pitaevskii equations (GPE) used for coherent 
ground-state of Bose condensed dilute atoms. This di¬ 
rectly leads to a wide range of analogies between exciton- 
polaritons and cold atoms. Actually, in semiconductor 
microcavity systems, we can investigate a wide range of 
physics including Bose-Einstein condensation mm and 
the collective quantum fluid nature [HHin]. On the other 
hand, one of the important properties of a semiconduc¬ 
tor system is the existence of dephasing [IIHIS], which 
induces decoherence. The investigation of the effect of 
decoherence, induced by excitonic dephasing, on the po- 
lariton dynamics is important both for understanding the 
physics of exciton-photon strongly coupled systems and 
for designing semiconductor microcavity devices such as 
nonlinear optical devices and polariton-based qubits. 

In this letter, we show that excitation induced dephas¬ 
ing (EID) plays an important role in the dynamics of 
polaritons in a semiconductor microcavity. The inves¬ 
tigation is experimentally carried out by time-resolved 
femtosecond pump-probe optical spectroscopy. Eor the 
theoretical description of our results, we utilize the ex¬ 
citonic Bloch equations (EBE) approach, taking into ac¬ 
count separately the coherent part of the polariton pop¬ 
ulation and an incoherent population of excitons [Hill]. 
We study the role of exciton-exciton interactions, photon- 


assisted exchange scattering, and EID effects on the po¬ 
lariton dynamics. The experimental results are very well 
reproduced by EBE and not by the exciton-photon GPE, 
which assumes that excitons are in a coherent limit. 

The experiment is performed in a high quality GaAs- 
based microcavity m at the cryogenic temperature of 
4K. A single 8 nm Ino. 04 Gao. 96 As quantum well is em¬ 
bedded in between two GaAs/AlAs distributed Bragg- 
reflectors. The Rabi splitting energy is 2D=3.45 meV at 
zero cavity detuning [9]. Eor the accurate measurement 
of the transmitted probe beam, we employ a heterodyne 
pump-probe setup with a degenerated beam configura¬ 
tion at /c = 0 /im' which dramatically increases 

the signal-to-noise ratio. The pump and probe pulses 
originate both from a broadband few hundreds femtosec¬ 
ond Ti:Sapphire laser. The center of the laser spectrum 
is set between the lower and upper-polariton peaks. Ad¬ 
ditionally, noise coming from laser spectrum envelope is 
removed with the aid of a numerical low-pass filter. The 
experimental setup is described in detail in our previous 
papers uni HI]. In order to avoid the complex effects 
of biexcitons HaHH], the pump and probe beams are co- 
circularly polarized. We obtain a time delay map through 
successive measurements of the pump-probe spectrum. 

The experimental results are presented in Eig. (a) 
and (c), showing the probe spectra as a function of 
pump-probe time delay. The cavity detuning is set at 
€c — €x = O.S meV, where ec(x) is cavity mode (exciton) 
energy. In this figure, the pump pulse arrives before (af¬ 
ter) the probe pulse at positive (negative) pump-probe 
delays. Eor low pump intensity (Eig. [^(a)), we observe 
two polariton branches (lower and upper) at both posi¬ 
tive and negative delays and the lower polariton shows a 
maximum blue-shift at zero delay. The delay dependence 
of the lower-polariton blue shift is asymmetric with re¬ 
spect to zero delay. While the blue shift gradually de¬ 
creases at negative delays, it stays at positive delays for 
long time. No clear energy shift of the upper polariton 
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resonance is seen for all delays. At a high pump intensity 
(Fig. [I] (c)), a triple peak structure appears at negative 
delays, while a single peak exists at positive delays. With 
the aid of numerical simulations based on EBE, we show 
that such behaviours originate from a long-living incoher¬ 
ent population and a short-living coherent polarization of 
excitons. 

Eor the analysis of the experiment, we use EBE [MlIIS]. 
The starting point of EBE is a bosonic exciton-photon 
Hamiltonian: 




H = J dx 

+n{ipl'tp^ + 'tPli’c) 


2mc 


+ j dx j dx' ^'tpllpl'Vey,{x-x')'tp:, 1 p'^ (1) 
-Vpae(x - + i’li’l'i’x'tpc) ■ 


ipx{c) arid are exciton (photon) field creation and 

annihilation operators. They obey the Bose commuta¬ 
tion relation, ['0cc(c)5'0^(c)] = — x'). This Hamilto¬ 

nian can be obtained from the electron-hole Hamiltonian 
via a boson mapping method called the Usui transfor¬ 
mation m- Since the exciton mass is large, the kinetic 
term of the exciton is neglected. The interactions are as¬ 
sumed to be contact interactions: V^x(x—x') = gS(x—x') 
and Upae(x — x') = gpae^(x — x') [191 ED]- The exciton- 
exciton interaction potential Wx is associated with the 
Coulomb exchange scattering. The term Upae is a pho¬ 
ton assisted exchange scattering [21] and contributes to 
the reduction of the Rabi coupling, which is the remi¬ 
niscence of the fermionic nature of the exciton [14]. In 
order to obtain a closed set of equations, we truncate the 
hierarchy by applying the following assumptions such as 

ii’x'ipx) = 0, and ~ =0. We 

define the exciton population as N{x^t) = and 

the exciton polarization as P{x,t) = Assuming 

factorization between the photon and exciton, we define 
E(x^t) = (V^c)- With the aid of the Heisenberg equation 
of motion, the EBE then reads El 

ihN = -iT^N - 2i(fi - 2gp^eN)lm[PE*] 
ihP = {ex + goN — ijx{P^)) P + (U — 2^paeA^)E (2) 
h‘^ 

ihE = (ec - -—- i'yc)E + (U - gps.eN)P - /ext- 
Zrric 

To obtain the above equations, the interaction constant g 
is phenomenologically divided into a real and imaginary 
part: g = go — ig'. The real part go is associated with 
an energy renormalization, while the imaginary part g' 
represents the strength of EID, which is also referred to 
as collisional broadening. The ratio of the constants is 
estimated as gpae/go — hCl/dUsEijao [22]. is the sat¬ 
uration density of excitons. and ao are respectively 
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FIG. 1. (color online) Measured and simulated probe trans¬ 
mission are shown as a function of energy and time delay 
between pump and probe pulse. The spectra are measured 
for two different pump intensities: 1.48x10^^ (1 mW) (a) and 
7.4x10^^ (5 mW) (c) photons pulse“^ cm“^. Corresponding 
simulated spectra are attached below the measured spectra 
((b) and (c)). The black dashed lines represent the lower 
and upper-polariton peak energies without pump pulse. The 
white dashed lies are the cavity photon and exciton energies. 
In the simulation, the intensity Iq is defined as Iq = 0.8/^o- 

the exciton binding energy and Bohr radius. The con¬ 
stants, jx{N) and Tx are respectively polarization de¬ 
phasing and population decay rate of excitons. In gen¬ 
eral, ^x{N) is written as [TTl [23l[24] . 

j,{N) = T,/2 + j:+g'N, 


(3) 
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where 7 * is the pure dephasing term. In the terminol¬ 
ogy of two level systems, and 7 ^, correspond to the 
inverse of Ti and T 2 times respectively . The EID con¬ 
stant g' introduces a phenomenological linear increase of 
the dephasing that depends on the exciton population N^ 
which plays an important role in our experiment. 

The advantage of EBE, compared to GPE, is that we 
can apply independent decay rates for the coherent polar¬ 
ization and incoherent population and calculate the time 
evolution of each. Indeed, in the commonly used GPE, a 
factorization, = (' 0 |)(' 0 a,) is implicitly assumed 

[20] . Therefore, the dynamics of the population A^(x, t) is 
uniquely determined by the polarization P(x, t) through 
N = \P\‘^. In the EBE, this condition is satisfied when 
neither pure dephasing nor EID exist: 7 * = 0 meV and 
g' = 0 meV (coherent limit). The EBE are close analogue 
of the optical Bloch equations (OBE) [l3l|25], however 
differing since OBE are based on a two-level electron-hole 
system, while EBE are based on a bosonic exciton basis 

m- 

To reproduce the experiments, Tx and 7 c are chosen 
to be 0.01 meV and 0.1 meV respectively. The pure de¬ 
phasing is set to 7 * = 0.1 meV [26], additionally, we 
include EID as the primary decoherence mechanism in 
our simulations. We set the interaction constants as 
g' = 0.4^0 and ^pae = O.S^fo* /ext is the excitation 
photon field and is assumed to be a Gaussian pulse: 

/ext = - to)^/('r^))exp(-iw7''''0- We 

set at the center of both polariton branches and 

set r=0.5 ps. 

Eor the calculation of the pump-probe dynamics, we 
use a coupled-mode approximation: A^(x, t) = + 

^ population is a real value), 

P(x,t) = + ppr^ik.^ _|_ pid^-ik.^^ P(x,t) = 

pp^ + ^ pid^-ik.^^ Pqp example, P^^, P^G and 

pid represent pump, probe, and idler component of the 
polarization respectively. Gonsidering the conservation 
of momentum, we obtain 8 coupled equations (See Sup¬ 
plementary material). The pump and probe pulses are 
introduced as and respectively and the trans¬ 
mitted probe signal is obtained through E^^. This is 
the standard method of calculating a transient four-wave 
mixing signal in optical Bloch equations EZIIIH]. Since 
the wave number of the probe is sufficiently small, we 
neglect the momentum dispersion of the photon mode. 

The simulated probe transmission spectra are given in 
Eig. (b) and (d) for two different pump intensities. 
There are striking similarities between the experimental 
and simulated spectra. Eirstly, Eig. [^(b) features a long 
lasting (^ h/Tx) blue-shift of the lower-polariton in the 
positive delay, while the blue-shift builds up on a shorter 
time scale (^ 2fi/(7c + 7 *)) in the negative delay. On the 
other hand, the energy shift of the upper-polariton is al¬ 
most zero because of the cancellation of the blue and red- 
shift contributions induced by go and ^pae respectively. 
The high-density simulation (Eig. [^(d)) reproduces both 
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FIG. 2 . (color online) Simulated probe transmission as a func¬ 
tion of energy and pump-probe time delay without EID or 
pure dephasing {g' = 7 * = 0 meV). The other parameters are 
the same as those used in the simulation of Fig. 1. Fig. 2 (a) 
and (b) respectively correspond to 1 mW and 5 mW pump 
intensities. 

the occurrence of three peaks at negative delays and of 
a single peak at positive delays. At negative delays, the 
signal shows a three peak structure, this is because the 
signal is temporally convoluted due to the finite lifetime 
of the probe pulse in the sample. Therefore the side peaks 
arise from the portion of the probe which is transmitted 
before the arrival of the pump pulse; however, the mid¬ 
dle peak builds with the arrival of the pump pulse in the 
short lifetime of the polarization given by ^ 2h/ ( 7 c+ 7 ^^)- 
At positive delays, since the probe polarization is always 
affected by the long lasting pump population, the single 
middle peak remains for a long time (^ h/Tx). If the 
pump intensity is further increased, eventually the term 
jx (N) becomes comparable to the effective Rabi coupling 
ft — 2 ^gpaeA^ and the central peak asymptotically reaches 
the cavity mode Cc, which is the signature of a strong to 
weak coupling transition [29]. 

For the purpose of better understanding the effect of 
the incoherent exciton population and of EID on the two 
polariton resonances, we present in Eig. [^a simulation 
without EID or pure dephasing {g' = 7 * = 0 meV). The 
other parameters are same as for Eig. The probe 
transmission of Eig. is a simulation of the coherent 
limit, where the polarization decay (dephasing) rate is a 
half of the population decay rate {^x{N) =fx/‘^ = 0.005 
meV). We find that the polariton branch is broadened to¬ 
wards the high energy side because of “dynamical energy- 
shift”. Namely, the mean-field energy-shift of polariton 
temporarily decreases following the decay of the polari- 
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FIG. 3. (color online) Simulated time evolution of the k = 
0 /im“^ polarization |P|^ and population N as a, function 
of time after an arrival of a single pulse. Simulations with 
EID and pure dephasing (a,c) and in a coherent limit {g' = 
7 * = 0 meV) (b,d) are presented for two different laser pulse 
intensities 1 mW (a,b) and 5 mW (c,d). The parameters 
are the same as the simulation of Fig. 1. The dashed lines 
represent scaled laser pulses. 


ton density. The time integration of the temporal de¬ 
crease of the energy-shift introduces a broadening of the 
polariton branches. The blue shift of both polariton res¬ 
onances decay two times faster at positive delays than 
they emerge at negative ones. In this limit, the dynamics 
of the exciton population is uniquely determined by the 
polarization and we can replace the population A/'(x, t) 
with the square of the polarization |P(x,t)p. Here the 
three sets of equations can be reduced to two equations 
composed of the exciton polarization and electric field, 
this is the commonly used exciton-photon GPE [201 [30]. 
Clearly, the GPE cannot reproduce the dynamics of po- 
laritons in the presence of EID, which is evidenced by 
the huge differences between Fig. [^and the experiments 
(Fig. 0 (a) and (c)). In particular, the high pump in¬ 
tensity simulation (Fig. (b)) reproduces neither the 
three peak structure nor the disappearance of the quan¬ 
tum beat pattern for the positive delay. This simulation 
also implies that the bleaching of the upper and lower 
polariton resonances in the positive delay (Fig. 0 (c)) is 
associated with EID. The transition to the weak coupling 
regime can also be observed with the GPEs, however the 
very strong dynamical blue-shift effects completely differ 
from the observed experimental behaviour. 

In Fig. we present simulated real time evolutions of 
exciton polarization P{t) and population N{t) at k = 0 
after the arrival of a single laser pulse. Without 
EID (coherent limit), the time evolution of the exciton 
population N{t) is found to coincide with that of |P(t)p, 
which supports the relation N{t) = \P{t)\‘^ and the fac¬ 
torization (iplipx) = ('0|)('0a,), which is assumed in the 
exciton-photon GPE. In this case, the lifetime of the sys¬ 
tem is mainly determined by the short photon lifetime 


(^ hj^c)- On the other hand, with EID, we have to dis¬ 
tinguish between the dynamics of the polarization and 
that of the population. While the exciton polarization 
P{t) is directly coupled to the cavity photon E{t) (not 
shown), there is no direct coupling between the exciton 
population N{t) and the photon E{t). Therefore, while 
the polarization decays with a lifetime of the same order 
as that of the cavity photon (^ the population de¬ 

cays independently and stays for a long time (^ h/Tx), 
even after the disappearance of the polarization. Intu¬ 
itively, the EID process can be understood as follows; 
microscopically, the exciton-exciton collisions introduce 
an energy fluctuation, which gives an additional random 
phase to the time evolution of the exciton field opera¬ 
tor '0a,(x, 0)6*^^+^^^)^/^. Since the phase Se is random, 
the expectation value P(x, t) = ('0a,(x,t)) shows a decay 
[3T] , which is the origin of the imaginary part of the inter¬ 
action constant g' (EID) [22 [32]. Meanwhile, the the en¬ 
ergy fluctuations affect neither the term '0|(x, t)'0a,(x, t) 
nor its expectation value A^(x, t) = ('0j(x, t)'0a,(x, t)) 
due to a phase cancellation m- Finally, we comment 
that the incoherent exciton population N should be in¬ 
terpreted as the “inactive” excitonic reservoir already dis¬ 
cussed in the context of non-resonantly excited polariton 
condensates [33l |34| . Actually, both incoherent exciton 
population featured in Fig. [^and inactive excitonic reser¬ 
voir have a long lifetime and nevertheless contribute to 
the energy shift of the polariton resonances [35] . 

In conclusion, we investigated the coherent dynamics 
of exciton-polaritons by femtosecond pump-probe spec¬ 
troscopy in a high quality semiconductor microcavity. 
We demonstrated that excitation-induced dephasing, as 
a manifestation of exciton-exciton interactions, is a ma¬ 
jor mechanism for the breaking of the strong coupling 
regime. We showed that the experimental results could 
only be simulated with the inclusion of EID in the exci¬ 
tonic Bloch equations. 

The present work is supported by the Swiss National 
Science Foundation under project N° 135003 and the Eu¬ 
ropean Research Council under project Polaritonics con¬ 
tract N° 291120. The polatom network is also acknowl¬ 
edged. 


E-mail: naotomo.takemura@epfl.ch 

[1] G. Weisbuch, M. Nishioka, A. Ishikawa, and Y. Arakawa, 
Phys. Rev. Lett. 69, 3314 (1992). 

[2] R. Gerna, Y. Leger, T. K. Parai'so, M. Wouters, 
F. Morier-Genoud, M. T. Portella-Oberli, and B. De- 
veaud, Nature comm. 4 (2013). 

[3] A. Amo, T. Liew, G. Adrados, R. Houdre, E. Giacobino, 
A. Kavokin, and A. Bramati, Nature Phot. 4, 361 (2010). 

[4] S. S. Demirchyan, I. Y. Ghestnov, A. P. Alodjants, M. M. 
Glazov, and A. V. Kavokin, Phys. Rev. Lett. 112, 
196403 (2014) 






















5 


[5] L. Dominici, D. Colas, S. Donati, J. P. Restrepo Cuar- 
tas, M. De Giorgi, D. Ballarini, G. Guirales, J. C. Lopez 
Carreno, A. Bramati, G. Gigli, E. del Valle, F. P. Laussy, 
and D. Sanvitto, Phys. Rev. Lett. 113, 226401 (2014) 

[ 6 ] J. Kasprzak, M. Richard, S. Kundermann, A. Baas, 
P. Jeambrun, J. Keeling, F. Marchetti, M. H. Szyman- 
ska, R. Andre, J. L. Staehli, V. Savona, P. B. Littlewood, 
B. Deveaud, and L. S. Dang, Nature 443, 409 (2006). 

[7] H. Deng, G. Weihs, C. Santori, J. Bloch, and Y. Ya¬ 
mamoto, Science 298, 199 ( 2002 ). 

[ 8 ] A. Amo, J. Lefrere, S. Pigeon, C. Adrados, C. Ciuti, 
1. Carusotto, R. Houdre, E. Giacobino, and A. Bramati, 
Nature Phys. 5, 805 (2009). 

[9] V. Kohnle, Y. Leger, M. Wouters, M. Richard, M. T. 
Portella-Oberli, and B. Deveaud-Pledran, Phys. Rev. 
Lett. 106, 255302 (2011) 

[10] V. Kohnle, Y. Leger, M. Wouters, M. Richard, M. T. 
Portella-Oberli, and B. Deveaud, Phys. Rev. B 86 , 
064508 (2012) 

[11] J. Shah, Ultrafast spectroscopy of semiconductors and 
semiconductor nanostructures, Vol. 115 (Springer Science 
& Business Media, 1999). 

[12] M. Kira and S. Koch, Progress in quantum electronics 
30, 155 (2006). 

[13] H. Wang, H. Q. Hon, and B. E. Hammons, Phys. Rev. 
Lett. 81, 3255 (1998) 

[14] G. Rochat, C. Ciuti, V. Savona, C. Piermarocchi, 
A. Quattropani, and P. Schwendimann, Phys. Rev. B 
61, 13856 (2000) 

[15] D. Wang and M.M. Dignam, Phys. Rev. B 79, 165320 
(2009) 

[16] R. Stanley, R. Houdre, U. Oesterle, M. Gailhanou, and 
M. Ilegems, Appl. Phys. Lett. 65, 1883 (1994). 

[17] N. Takemura, S. Trebaol, M. Wouters, M. T. Portella- 
Oberli, and B. Deveaud, Phys. Rev. B 90, 195307 (2014), 

[18] N. Takemura, S. Trebaol, M. Wouters, M. T. Portella- 
Oberli, and B. Deveaud, Nature Physics 10, 500 (2014). 

[19] L. Pitaevskii and S. Stringari, “Bose-einstein condensa¬ 
tion (international series of monographs on physics),” 
(2003). 

[20] 1. Carusotto and C. Ciuti, Reviews of Modern Physics 
85, 299 (2013). 

[21] M. Combescot, M. Dupertuis, and O. Betbeder-Matibet, 
EPL (Europhysics Letters) 79, 17001 (2007). 

[22] D. Sarchi and V. Savona, Phys. Rev. B 75, 115326 (2007). 

[23] T. Baars, M. Bayer, A. Forchel, F. Schafer, and J. P. 
Reithmaier, Phys. Rev. B 61, R2409 (2000) 

[24] C. Ciuti, V. Savona, C. Piermarocchi, A. Quattropani, 
and P. Schwendimann, Phys. Rev. B 58, 7926 (1998) 

[25] X. Li, T. Zhang, C. N. Borca, and S. T. Cundiff, Physical 
review letters 96, 057406 (2006). 

[26] A. Honold, L. Schultheis, J. Kuhl, and C. W. Tu, Phys. 
Rev. B 40, 6442 (1989) 

[27] T. Yajima and Y. Taira, Journal of the Physical Society 
of Japan 47, 1620 (1979). 

[28] D. B. Turner, P. Wen, D. H. Arias, and K. A. Nelson, 
Physical Review B 84, 165321 (2011). 

[29] F. Quochi, G. Bongiovanni, A. Mura, J. L. Staehli, B. De¬ 
veaud, R. P. Stanley, U. Oesterle, and R. Houdre, Phys. 
Rev. Lett. 80, 4733 (1998) 

[30] Strictly speaking, the coupled-mode approximated EBE 
in coherent limit is slightly different from a coupled-mode 
GPE. 

[31] S. Schmitt-Rink, S. Mukamel, K. Leo, J. Shah, and D. S. 


Chemla, Physical Review A 44, 2124 (1991). 

[32] C. Cohen-Tannoudji, J. Dupont-Roc, G. Grynberg, and 
P. Thickstun, Atom-photon interactions: basic processes 
and applications (Wiley Online Library, 1992). 

[33] K. G. Lagoudakis, F. Manni, B. Pietka, M. Wouters, 
T. C. H. Liew, V. Savona, A. V. Kavokin, R. Andre, 
and B. Deveaud-Pledran, Phys. Rev. Lett. 106, 115301 
( 2011 ) 

[34] M. De Giorgi, D. Ballarini, P. Cazzato, G. Deligeorgis, 
S. 1. Tsintzos, Z. Hatzopoulos, P. G. Savvidis, G. Gigli, 
F. P. Laussy, and D. Sanvitto, Phys. Rev. Lett. 112, 
113602 (2014) 

[35] L. Perrier, E. Wertz, R. Johne, D. D. Solnyshkov, 
P. Senellart, 1. Sagues, A. Lemaitre, G. Malpuech, and 
J. Bloch, Phys. Rev. Lett. 106, 126401 (2011) 


SUPPLEMENTAL MATERIAL: COUPLED MODE 
EQUATIONS FOR PUMP-PROBE DYNAMICS 

We present explicit forms of the coupled mode equa¬ 
tions for simulating pump-probe spectra. We restrict in¬ 
volving modes into three modes: pump, probe and idler. 
The polarization P, population N and electric field E are 
respectively written as N = 

p = + pprf,ikr pid^-ikr^ p ^ ppu ppr^ikr 

pid^-ikr ^ Substituting these representations into the ex- 
citonic bloch equations (EBE) in the manuscript and ne¬ 
glecting components such as and we obtain 

8 coupled equations of motions. Eirstly, the population 
N read, 

ihNP'^ = -iV^NP'^ 

+ 2 flpae(ap„VP®® + buNP^ + bprNP^*) 

ihNP^ = -ir^NP^ 

-nbpr + 2gpUbprNP^ + apuNP^ + cNP^*) 


Now the quantities ap„, bpr, bid and c are given by 


apu = 2 iIm(PP®‘£’P’®* -h pP^EP^* -|- p^dpid*) 


bpr = PP^E 


\id^ 


^ — ppr -^id^ 


pid^ ^pu _|_ pprppu^ ppu^ ppr 

ppr^ ppu _j_ pidppu^ _ ppu^ pid 

pid^ ppr 


The equations of motion of exciton polarization is written 
as 

ihPP^ = {e^ - i'y^)PP^ + giNP^PP^ + NP^*PP^ + NP^P^<^) 

- 25pae(VP’®£’P’® -h NP^*EP'^ + NP^'E^'^) 

ihPP^ = - z 7 a,)P^’' -h g{NP^PP^ + jqprppu^ 

-hCPP®' - 25pae(VP’®PP®' -h NP'^EP^) 

= (e^ - i7a;)P*‘^ -f g{NP^P^'^ + NP^*PP^) 

- 25pae(VP®‘P®‘' -h NP^*EP'^). 
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Finally, the electric fields follow the following equations: 
ihEP^ = (ee - 

+ NP^P^'^ + NP^*PP^) - /I*" 

ihEP^ = (ee - hc)EP^ + 

-gp.eiNP^P^’' + MP’-pp-) - fi:, 

ihE^'i = (e^ - Z7e)-E’*‘' + - 5pae(A^^“P*‘' + NP^*PP^). 



